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Viscous Corrections to Spectra, Elliptic Flow, and HBT Radii 
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^ , I compute the first viscous correction to the thermal distribution function. With this 
correction, I calculate the effect of viscosity on spectra, elliptic flow, and HBT radii. 
Indicating the breakdown of hydrodynamics, viscous corrections become of order one for 
Pt ~ 1.5 GeV. Viscous corrections to HBT radii are particularly large and reduce the 
^ . outward and longitudinal radii. This reduction is a direct consequence of the reduction 
! in longitudinal pressure. 
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O ' Ideal hydrodynamics describes a wide variety of data from heavy ion collisions [ 0, • 
In particular ideal hydrodynamics successfully predicted the observed elliptic flow and its 
dependence on mass, centrality, beam energy, and transverse momentum. Nevertheless, 
O I the hydrodynamic approach failed in several respects. First, above a transverse momen- 
^ I tum pr ~ 1-5 GeV the particle spectra deviate from hydrodynamics and approach a power 
law. Second, HBT radii are significantly too large compared to ideal hydrodynamics [ H . 
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1. Viscous Corrections 



Considering the partial success of ideal hydrodynamics, viscous corrections may provide 
a natural explanation for these failures. 
5^ i For an ideal Bjorken expansion, the entropy per unit space-time rapidity (rs) is con- 
served. For a viscous Bjorken expansion the entropy per unit rapidity increases as a 
function of proper time [ ffl 
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where r] is the shear viscosity. In this equation and below we have neglected the bulk 
viscosity. For hydrodynamics to be valid, the entropy produced over the time scale of the 
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expansion (to wit, t^) must be small compared to the the total entropy (rs). This leads 
to the requirement that 

- « 1 , (2) 
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where we have defined the sound attenuation length as F^ = Perturbative estimates 
of the shear viscosity in the plasma give ^ ~ 1. Below we take ^ = |, assuming that 
non-perturbative effects shorten equilibration times. 

Viscosity modifies the thermal distribution function. This modification influences the 
observed particle spectrum and HBT correlations. The formal procedure for determining 
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the viscous correction to the thermal distribution function is described in references [ ^ . 
However, the basic form of the viscous correction can be intuited without calculation. 
First write f{p) = /o(l + g{p)), where fo{^) is the equilibrium thermal distribution 
function and g{p) is the first viscous correction. g{p) is linearly proportional to the 
spatial gradients in the system. Spatial gradients which have no time derivatives in the 
rest frame and are therefore formed with the differential operator = [Qf^i, — u^u,^)d'^ . 
For a baryon free fluid, these gradients are Vq,T, Vq?/", and {VaUp), where (VaUp) = 
+ VjsUa — | A^/? V-y?/'*'. Vq,T cau be eliminated in favor of the other two spatial 
gradients using the condition that T^^^Uy = eu^ and the ideal equations of motion. Vo-ti" 
leads ultimately to a bulk viscosity and will be neglected in what follows. Finally, (VaUp) 
leads to a shear viscosity. If g{j)) is restricted to be a polynomial of degree less than two, 
then the functional form of the viscous correction is completely determined, 

/ = /o(l + ^pV(V„n^)). (3) 

For a Boltzmann gas this is the form of the viscous correction adopted in this work. For 
Bose and Fermi gasses the ideal distribution function in Eq. ^ is replaced with foil ± /o) 

[i- 

The coefficient C is directly related to the sound attenuation length. Indeed, using the 
distribution function in Eq. ^ to determine the stress energy tensor, yields a relationship 
between the shear viscosity r] and the coefficient C, 

^ r^,U ^ ^ = jdV-^ foil + i^aUp)) . (4) 

For a Boltzmann gas, Eq. ^yields C = ^. 

The thermal distribution function is now completely determined. In the next section 
this correction is used to calculate corrections to the observables used in heavy ion colli- 
sions. 

2. Corrections to Spectra, Elliptic Flovi^, HBT Radii 

To quantify the effect of viscous corrections on spectra and HBT radii, I generalize the 
blast wave model. In the blast wave model used here, the matter undergoes a boost in- 
variant Bjorken expansion and decouples at a proper time of Tq = 6.5 Fm at a temperature 
of To = 160 MeV. The matter is distributed uniformly up to a radius of Ro = 10.0 Fm, 
with velocity profile up to a maximum velocity Uo = 0.5 cQ. The blast wave model with 
these parameters closely models the output of a full hydrodynamic simulation [ [l|, || and 
gives a reasonable fit of the data. 

The spectrum of produced particles is given by the Cooper-Frye formula, 

dNo + 6dN = j p^dT.^fo + 6f . (6) 

Fig. [l| shows the ratio of the correction compared to the ideal spectrum, 
^In contrast to common practice w'' = ^v^ is linearly rising: — 
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Figure 1. Viscousity corrected spec- Figure 2. Elliptic flow as a function trans- 

trum relative to uncorrected spectrum. verse momentum. The blast wave param- 

The band indicates where hydrodynamics eters (see text) are chosen to approximate 

breaks down. a AuAu collision at & = 6 Fm. 



To understand this figure qualitatively, consider a Bjorken expansion of infinitely large 
nuclei. The longitudinal pressure is reduced pl = p — Because the shear tensor is 
traceless, the transverse pressure is increased, Pr = P + Thus, the matter distribution 
is pushed out to larger pt by the shear in the longitudinal direction. More mathematically, 
the ratio of the corrected spectrum to the uncorrected spectrum is given by, 

dNo At\\t) \T ) 2 [Ki{^) ) 

For large px we find, ~ ^ ' reproduces the shape and dependence of the 
full viscous blast wave calculation shown in Fig. |I|. 

Viscous corrections become of order one when the pt of the particle approaches 1.4 GeV. 
This signals the breakdown of the hydrodynamic approach. In fact, ideal hydrodynamics 
generally fails to reproduce the single particle spectra above pt of 1.5 GeV. Viscosity 
provides a ready explanation for this breakdown. 

In non-central collisions elliptic fiow is calculated using the spectrum indicated in Eq. 
^ In non-central collisions, the matter is assumed to have a cylindrical distribution but 
the fiow velocity has an elliptic component, u'^{r,(j)) = ul^{l + 2u2Cos{(p)). For non- 
central collisions the parameters are: Uo = 0.5, U2 = 0.1, Ro = 6 Fm and Tq = 4.0 Fm. As 
illustrated in Fig. viscosity reduces elliptic fiow by a factor of three. 

Next consider viscous corrections to HBT radii. The HBT radii are calculated with the 
method of variances. First the ideal radii parameters are displayed in Fig. ^. The results 
are typical of the blast wave parametrization. Next the viscous correction to the blast 
wave results are illustrated in Fig. H. 
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Figure 3. HBT radii for the ideal blast 
wave model described in the text. 
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Figure 4. Viscous corrections to the ideal 
HBT radii illustrated in Fig. ^. 



Viscous corrections to and Rq are large and negative. This may be understood 
qualitatively by again considering a simple Bjorken expansion of infinite nuclei. The 
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longitudinal pressure is reduced, Pl = P — Therefore the Pz distribution 
narrower. However, by boost invariance the single particle distribution (^i^) is a function 
of \y — r]\, which yields the relation, 

dN 
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dydz 
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The z distribution (^^) at mid rapidity is therefore narrower because the longitudinal 
pressure is reduced. This implies that -R| = (z^) decreases due to the viscous longitudinal 
expansion. In summary, viscosity provides a simple explanation for the very large radii 
predicted by ideal hydrodynamics. 
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